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We show that a spherically symmetric perturbation of a dust dominated £1=1 FRW universe in 
the Newtonian gauge can lead to an apparent acceleration of standard candles and provide a fit 
to the magnitude-redshift relation inferred from the supernovae data, while the perturbation in 
the gravitational potential remains small at all scales. We also demonstrate that the supernovae 
data does not necessarily imply the presence of some additional non-perturbative contribution by 
showing that any Lemaitre-Tolman-Bondi model fitting the supernovae data (with appropriate initial 
conditions) will be equivalent to a perturbed FRW spacetime along the past light cone. 
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INTRODUCTION 



A basic principle which has enabled substantial progress in cosmology, referred to as the Copernican Principle, is 
the assumption that no observer has a privileged position in space. Coupled with the observed isotropy of the universe 
around us, this principle leads to the conclusion that the universe is homogeneous. Although the consistency of the 
current cosmological models strongly supports this idea, the Copernican Principle has never been explicitly tested; 
,Ih ■ even the possibility of such testing has only been raised rather recently [l|, S 3 • Even if the Copernican Principle 
Oh- is in general valid, there is still the possibility that it is mildly violated at some cosmologically significant level and 
that observers could occupy inequivalent spatial locations. In that case, their observations might mislead them into 
H ' inferences about the universe as a whole and its global parameters which are incorrect. 

Doubts about the validity of the Copernican Principle lie very much at the heart of some of the explanations 
1^1 ' put forward to account for the apparent acceleration implied by the Type I supernovae data 0]. As is well known, 
all cosmological data is consistent with a homogeneous ACDM concordance model, which assumes the existence of 
a cosmological constant A that makes up about 70% of the current energy budget of the universe. Equally well 
acknowledged is the fact that, from a theoretical point of view, the incredibly small value of the required cosmological 
constant, about 10 -120 M p i, is a mystery. Moreover, there is also a peculiar coincidence that A appears to have 
become important at the time when structures started to form. This may be taken as a hint towards the possibility 
that inhomogeneities could in some fashion be responsible for the apparent acceleration. One school of thought 
^ — I , attempts to explain it as the influence, sometimes called backreaction, of averaging inhomogeneities on the average 
expansion of the universe [H, @, 0, S, @, E3, U , E EH • Backreaction would be an inherently nonlinear effect due 
O s \ to the nonlincarity of general relativity. In this case the Copernican Principle could still hold in some statistical 
sense. Another possibility is to abandon the Copernican Principle altogether and assume that our location in the 
^ , universe is exceptional, thereby fooling us into thinking that a cosmological constant exists, while in reality it is only 
a trick of the light. Such considerations have mostly been put forward within the context of spherically symmetric 
Lemaitre-Tolman-Bondi (LTB) models 14, 15. Til 17 [. Many authors have considered various possibilities offered 
by LTB models [Hj], which include inhomogeneities in the matter distribution as well as inhomogeneities in the 
expansion rate; an example is a dust universe with a local underdensity [3] that seems to be able to reproduce the 
SNIa observations. 

There is an ongoing debate about the role of inhomogeneities, particularly whether the effects of inhomogeneities are 
too large to be describable within the standard perturbed FRW description. In contrast to various model calculations, 
in the present paper we do not make an attempt at model building and full-blown data fitting. Rather, we address 
the generic question of whether the SNIa data requires a deviation from the Copernican Principle that is non- 
perturbative. To explore this issue, we consider an FRW spacetime that has small perturbations in the gravitational 
potential ip (although the density perturbations could well be large). We assume an SI = 1 dust universe and, 
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for simplicity, spherically symmetric perturbations, and solve for ip and the velocity field v such that the distance- 
luminosity relation of the SNIa data is satisfied; this we do in Newtonian gauge. We find that the SNIa data can 
be fitted with tp and v remaining small during the whole evolution, meaning that the perturbed FRW framework 
works well and the Copernican Principle is thus only weakly violated in the sense that there are only small deviations 
from homogeneity. Furthermore, one might then argue that any effects of averaging over these small inhomogeneous 
perturbations, i.e. the backreaction, remain negligible (20j ; due to the nonlinear nature of backreaction, the effects 
of averaging can only arise at second order in the perturbations and would hence be extremely small. Indeed, if the 
spacetime metric of the real universe were a perturbed FRW, it would imply that the backreaction is insignificant 
and thus could not solve the dark energy problem 1 3. l2ll However, we should note that the backreaction issue is 
somewhat more involved, and refer the reader to 0, 1 lOL l2l| and the references therein. Indeed, in a nonlinear theory 
there could arise additional non-perturbative contributions to metric perturbations which are not manifest at any 
order of perturbation theory performed about the FRW background but could be revealed by a careful comparison of 
data and results obtained from perturbation theory. As we will demonstrate, at least the supernovae data does not 
necessarily imply the presence of some additional non-perturbative contribution. 

Finally, we contrast our result with LTB models, which under certain conditions are known to be physically 
equivalent to perturbed FRW models 0, [1(1 HH j although it has also been claimed that the LTB models capable of 
accounting for the SNIa data cannot be obtained from a perturbed FRW model Here we show explicitly that 
such conditions arc fulfilled for realistic LTB models fitting the SNIa data and thus conclude that an explanation 
of the SNIa data in terms of (spherically symmetric) inhomogeneities does not require non-perturbative effects that 
would not be already present in a perturbed FRW description. 



A SPHERICAL PERTURBATION 



Let us assume that the spacetime can be described by a perturbed FRW model in the Newtonian gauge (vanishing 
anisotropic stress) 



ds 2 = a 2 [-(1 + 24>)drf + (1 - 2^)dx 2 ] 



(1) 



The effect of linear perturbations on the luminosity distance has been calculated in the past [23|, |2J, |25j . We follow 
here the notation of Bonvin, Durrer and Gasparini [25| . who find the luminosity distance di(r7 s ,n) of a source at 
conformal time ry s in a universe described by the metric |T]) as given by [25| 
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where n = — e r refers to the spatial direction of light propagating towards the observer, rj s is the conformal time 
at emission, and Ar) = r\ — rj s . The integrands of the above integrals are taken to depend on A (A) as well as on 
x = x G — n (rj — A), i.e. they are evaluated on the past light cone. The subscript "o" refers to the observer today. 

In Eq. @ rj s refers to the background conformal time corresponding to the emission of light from the source. 
However, in the presence of perturbations a fixed value of rj does not correspond to a spatially homogeneous value 
of z. Furthermore, it is z, not 77, which is measurable and by considering any quantity as a function of redshift in 
effect one slices the universe in slices of constant z, not 77. From Eqs. (56) and (57) of [25j we find for the luminosity 
distance as a function of redshift the expression: 
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from which we obtain 
d L {z,h) 



At? = [ 2Ai] + — ] (v • h- ij) ) - ( 3Ar/ + ] v s • n + [ A 77 + -i 



i + z ' v w»/ v ro/ \ ' nj \ ' n 

+2 J + 2 \2Ar]+ ^-J / dAn-V?/> 

+2 y dA y d\n • Vi/} - J d\ J dX(X- r] s ) (VV - n • V (n • W)) • (4) 

Eq. ([4]) can be used to determine the fluctuations of the luminosity distance induced by a Gaussian perturbation 
(25| , such as the one commonly assumed to arise from inflation and describing structure formation in the concordance 
cosmology. However, we will turn the reasoning the other way round and ask: Assuming the lhs of Eq. (|4| is a 
measured function from SNIa observations and that the underlying model for the background universe is flat and 
matter dominated with VLm = 1 and A = 0, can we determine the perturbation tp (or v) needed to give the observed 
luminosity distance? In other words, can a local gravitational perturbation fool us into thinking that light has 
propagated in a universe with the estimated value for the cosmological constant? We will answer this question by 
assuming that the perturbation ip exhibits spherical symmetry for simplicity. As discussed in the Introduction, we 
are interested in the case where the background can still be described by the FRW metric. Furthermore, assuming 
linear theory, one can relate the velocity field v with respect to the Einstein-de Sitter background and the Newtonian 
potential by 

v = -J^W>. (5) 

Our goal then is to determine v and tp, given the lhs of (J4j) from observations (The d^ — z relation corresponding to 
an Q,\ ~ 0.7, CIm ~ 0.3 universe). 

It is convenient to derive a second order equation for the radial velocity v r , by acting with two derivatives along the 
past light cone ^ = d v + n- V, assuming that it is the dominant component. Indeed, if the radial component is to be 
able to account for the observed dj, — z relation, it will have to be be much larger than the other velocity components 
related to the standard small Gaussian density fluctuations. Given v r the potential can be determined via Eq. ([5]). 

Let us now proceed in the manner just described. Taking two derivatives of (j4]) we obtain 
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where v sr = — v s • n. To derive this result we have used the following: For a matter dominated universe = \ and 
drjtp = 0. The latter relation allows us to write dip/drj = n • = — ■ v - To eliminate the V 2 ip term we assumed 
that v r is dominant, i.e. that the main component of the peculiar velocity is radial. Thus 



Furthermore, 



Since = 0, we have that 



from which we get (r = Aij) 
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We can now express Eq. ([6]) in terms of the redshift z using 

dz 



H = HqvT+z , and drj = — - 
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(l+z)v sr = 



2 r 



1 + z 



(11) 



(12) 



This is our main result and can be seen as the equation which defines the perturbation in the (radial) peculiar velocity 
given an observed dh — z relation, acting as a source on the rhs. Given appropriate boundary conditions, its solution 
determines the velocity needed to give rise to the observed dh — z curve. We stress again that the rhs is zero for 
a matter dominated FRW universe and hence any deviation from zero is attributed to a perturbation. Thus, the 
physically relevant boundary conditions for our case are v r (z = 0) = -j^v r (z = 0) = 0. Allowing for curvature in the 
background model would bring complications to the perturbation theory and thus would require a separate analysis, 
so for simplicity and also motivated by the observations we restrict ourselves to the flat model. 
Using Eq. (|11|) we can calculate the rhs of Eq. (|T2|) to find 
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Assuming that the observed dh is described by the corresponding theoretical relation for a ACDM universe 

1 + 2 
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where a ~ 0.3, (3 ~ 0.7 and # ACDM = 70 km/s/Mpc, we have 
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In Fig. [T] we plot the solution to the equation of motion for the velocity field, given by Eq. (p~2|) . that gives rise to the 
luminosity distance given in Eq. (|15[) . We also construct the gravitational potential, assuming that tp is zero at large 
redshifts. We see that a spherically symmetric linear perturbation can indeed explain the observed acceleration with 
small velocities and small gravitational potentials. Note that the value H in the hypothetical fl = 1 dust universe 
is not really an observable since the model is unlikely to fit all cosmological data. For the purposes of the present 
approach, it is an adjustable parameter. Here we have chosen Ho = 50 km/s/Mpc for illustrative purposes only. As is 
obvious from the gravitational potential of Fig. [I] physically the perturbation corresponds to a spherical underdensity 
around us. 

The density profile on the light cone can be immediately obtained via the poisson equation 



4ira 2 G5p = V 2 ?/> - 3W 2 ^ 

Using Eqs. (TIT))) and (jTTj) we thus have for the total matter density p = 
the light cone: 



(16) 

p + Sp, where p = 3H$/(8nG)(l + z) 3 , along 
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TRANSFORMATION TO LTB COORDINATES 



Let us now turn to LTB models. The LTB line-element with the spatial origin at the symmetry center reads as 

ds 2 = -dP + {A ' {t \ T) } 2 dr 2 + A 2 {I f) (d9 2 + sin 2 6d^ 2 ) , (18) 
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FIG. 1: The velocity field v ar (z) and the gravitional potential tp(z) of the perturbed FRW universe. The plots on the right are 
close-ups of those on the left for redshifts z < 10. 




FIG. 2: The density profile p(z)/p cr it(z) (the red line). The green line corresponds to the background Einstein-de Sitter 
p(z) I pciit(z) — 1. The plot on the right is a close-up of that on the left for redshifts z < 10. 



where A(i, f) is the scale function having both temporal and spatial dependence and k(f) is a function associated 
with the curvature of t = const, hypersurfaces. Both t and A have a direct physical meaning, the first one being 
the proper time of observers comoving with the dust and the second the physical radius of a sphere with area Air A 2 . 
The coordinate f is an arbitrary label, devoid of physical meaning. We use the following shorthand notations for the 
partial derivatives: ' = and ' = J^. The LTB metric ([T5]) is an exact solution of the Einstein equations and the 
perfectly homogeneous FRW model is a special case, obtained in the limit: A(t, f) — > a(t)f and k(f) — » kr 2 , where 
a(t) is the FRW scale factor and k is the curvature constant. 

It has been shown that under certain conditions, the LTB metric (jT5J can be brought to the perturbed FRW form 
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([T]) via the following nonlinear gauge transformation [l(| , j 

A(t,f) 



(!+£(*») (19) 



(l+^(f(z),f(z)))=2iJ ' 1 (l-l/ % /TT^) . (21) 



a(t) 

t = t + £°(f,f) , (20) 

where a(f) = (t/t ) 2 ^ 3 is a fictitious FRW scale factor corresponding to an Q = 1 matter dominated universe. If 
the LTB model is to be close the the FRW, the functions f(i,f) and £°(t,f) must satisfy < 1. We will 

demonstrate that these conditions are indeed satisfied for our case on our past light cone. 

Before proceeding we should note that Eq. (|20| provides a relation between the clock times of the LTB and 
(fictitious) FRW observers. Their difference can be intuitively understood in two ways: In a local description, the 
LTB coordinates are synchronous, meaning that they are defined by the world-lines of the dust particles which, 
since pressure is absent, are geodesies. Therefore, their clocks will tick faster than those of any other non-geodesic 
coordinate system such as the FRW coordinates. The same conclusion is reached in a global description where the 
LTB spacctimc describes a physical underdensity around our position, meaning that the LTB clocks will tick faster 
with respect to the FRW coordinates in this case. This effect will be more pronounced as we move to the centre 
and thus is expected to decrease with decreasing redshift. Eq. (|19p defines the FRW distance r in terms of the 

physical distance A, a prescribed FRW scale factor a(t) = (t/to) 2 ^ 3 and a function £ which is chosen to ensure that 
on the light cone r takes the standard FRW form in terms of redshift: 

a(z) 

Furthermore, we choose the arbitrary f such that, on the light cone: 

d f A(z) = a(z) . (22) 

Let us now determine and £. In the LTB model, the angular diameter distance is related to the metric simply 
by 

d A (z)^A(r(z),t(z)) = 7 ^L , (23) 

(l + zy 

where the last equality is due to Etherington's theorem 26|, 27| stating that for geodesic light in any spacetime the 
luminosity and angular diameter distances are related through d L = (1 + z) 2 d,A- If the perturbed FRW and LTB 
metrics are to describe the same spacetime, just written in different coordinates, the observable distance-redshift 
relations should be the same: d A ^ :B = d^ CDM . So, by having the observable distance-redshift relations we can readily 
read off A(t, f) along the past light cone of the LTB model, that is A(z). £ can then be determined directly from (|21[) . 
Comparing the angular parts of metrics ([1]) and (p~8|) we see that, to leading order: 

A(t, f) ~ a(t)r [1 - <p(t, r)] ~ a(t)r [l - f) + H$)?(t, r)] , (24) 

where we have used (f2Tj|. In the previous Section we determined the v(z) (or ijj(z)) profile giving the ACDM luminosity 
distance in the perturbed FRW dust universe. Using Eqs. p^|) . (|2~lj) and (|2"3"|) , we can thus obtain the small quantity 
H£° from Eq. along our line of sight: 

M*) - TT7T~, — \ f 1 - -7^=1 ( 1 ~ V>W + H (z)e(zj) (25) 



H (l + z)\ VTTz~ 
=► g(z)g°(z) = d f r f (1+ 1 z) f (26) 

where <1a is obtained from Eq. ([T4")) [d,A = (1 + z)~ 2 dL). 

The outcome is shown in Fig. [31 where we display the shift between t and t in units of the inverse Hubble rate as 
well as £. As can be seen both H£° and £ are less than one. In principle, we could have kept higher orders in our 
expansions, but other than changing the numerical values of by O(30)%, we do not expect any qualitative change 
in this picture. Thus, the conclusion is that an LTB model fitting the SNIa data can be mapped into a perturbed 
FRW model by a nonlinear coordinate transformation of the form (|19j) . (I20|) . at least on the light cone. 
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FIG. 3: The functions H£°(z) and £(z) relating the perturbed FRW and LTB clock times. The plots on the right are close-ups 
of those on the left for redshifts z < 10. 

In the literature, there has been discussion about the so called "inverse-problem", that is, finding a map from 
a given luminosity distance to the corresponding LTB model [28], |2^, [j}(J. In fact, in [3(| an LTB model, whose 
cLl — z relation is equivalent to that of the ACDM model was constructed. Not so surprisingly, the resulting LTB 
model was found to be such that there is a void whose symmetry center is at the observer's position - assuming the 
early universe was homogeneous. As expected, our results are consistent with this: In the previous Section we found 
that the perturbed FRW fitting the ACDM — z curve contains a void around our location, given the universe is 
asymptotically Einstein-de Sitter. Since the LTB and perturbed FRW are the same spacetime, usually just written 
in different coordinates, the LTB model must also represent this void. So although we have not explicitly calculated 
the quantitative properties of the resulting LTB model, the physical content is known from a qualitative point of view 
and consistent with [3(|, as it should. 

Before summarizing our results, a few comments on the paper [ll[ are in place, as their results seem to contradict 
ours. As a toy model, they consider an LTB model consistent with the SNIa data with a large (1 Gpc) void in a 
universe approaching spatially to Einstein-de Sitter. They find that such models cannot be described in terms of a 
conformal Newtonian metric perturbed about an FRW background during the entire evolution without large peculiar 
velocities. However, our result is that the SNIa data can indeed be explained by the perturbed FRW model as well 
as by an LTB model, with these two descriptions being equivalent and with the peculiar velocity remaining small. 
Correspon ding ly, we find that the shifts between FRW and LTB coordinates are also small, in agreement with the 
analysis of [22J (see also |9|). In their analysis these authors include nonlinear but leading-order terms in the velocities 
which are absent from [111 ]. Indeed, since the solution is in the nonlinear regime, one would have to employ the full 
nonlinear transformation. However, more important is the fact that the peculiar velocities of the model considered 
in [ll| become large at redshifts z > 500, when the void represents a nonlinear superhorizon density fluctuation, so 
the model is not even in principle describable within the weak field description of Newtonian gauge. This means 
that there actually is no discrepancy between their consideration and ours: one can achieve large deviations from the 
perturbed FRW picture, but only by choosing initial conditions which seem observationally unrealistic. Indeed, when 
confronted with the CMB observations, a universe dominated by a Gpc sized void seems unlikely. 
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DISCUSSION 



We have established that the FRW fl = 1 dust universe with spherically symmetric perturbations in the gravitational 
potential tp(z) can fit the SNIa data. We found that physically the perturbation corresponds to a void around us. We 
do not pretend to know what could be the origin of the perturbation or the associated velocity field, and in fact it is 
rather obvious that the model as such is likely to be way too simplistic to explain all the relevant cosmological data. 
However, the interesting point here is that the SNIa data can be explained by small deviations from FRW; both ip(z) 
and the velocity field v(z) remain perturbative at all distance scales. Moreover, such a perturbative approach seems to 
capture all the essential features of the data. As we demonstrated, if the non-perturbative, spherically symmetric LTB 
model fits the SNIa data, it can be mapped into the perturbed FRW model by a nonlinear coordinate transformation, 
at least on the light cone. 



Our result seem to contradict where the claim is that the LTB metric cannot be expressed as a perturbed 



conformal Newtonian metric during the entire evolution of the universe without the existence of a peculiar velocity 
larger than the speed of light. However, they use a 1 Gpc LTB void model as a tool to demonstrate this and indeed, 
such case indicates the breakdown of the perturbative description: their model approaches Einstein-de Sitter only 
spatially and there were strong density contrasts in the early universe. The question is then whether such large 
nonlinear inhomogeneities are realistic. Observations seem to suggest that the inhomogeneities at the last scattering 
surface do not lead to such large underdense regions. Our model approaches asymptotically to Einstein-de Sitter along 
the past light cone, which implies that the early universe was very homogeneous and as our calculations explicitly 
demonstrate, the explanation of the supernovae data does not require a breakdown of the weak field description. As 
shown in Fig. [U the peculiar velocity remains small during the entire evolution along the past light cone and the 
mapping from LTB to perturbed FRW is possible. This is consistent with the results of (jl. [22|. 

We should like to stress that in fact we have not specified the LTB model quantitatively. That is, we have neither 
determined the LTB curvature function k(f) nor the mass function F(f), nor the entire evolution of A(t, f) inside our 
past light cone. However, the underdense bubble around us, represented by the perturbed FRW model, has to be also 
the physical content of the resulting LTB model from a qualitative point of view since, as we verified, the LTB and 
perturbed FRW are the same spacetime on our past light cone. In principle, we should have all the information to 
unambiguously determine also the quantitative properties of that LTB model: Our physical situation is such that the 
universe approaches asymptotically Einstein-de Sitter along the past light cone, so from this consideration it is obvious 
that the LTB model contains only growing mode inhomogeneities, so that the early universe was homogeneous and 
the Big Bang was homogeneous or the age of the universe is independent of the spatial coordinates: to(f) = constant. 
This requirement sets a constraint between the two boundary condition functions of the LTB model, thus leaving only 
one free function and a free parameter to, just like in the perturbed FRW case - the potential ip(z) (or alternatively, 
the velocity field v(z)) and a free parameter Hq. However, our interest here was not to concentrate on the exact model 
building. In fact, the issue of finding an LTB model that fits to a given d^ — z curve has been previously addressed by 
a few authors [28|, [29 . [3ll | and indeed, the physical content of the LTB model described here is consistent with the 
result obtained in [3fJ, where an LTB model was constructed by demanding its luminosity distance to be equivalent 
with the one in the ACDM model - their result is a void, the center of which we occupy. We hope to return to the issue 
of quantitatively determining the LTB model from luminosity distance data in future work. It will also be interesting 
to see whether the perturbative nature of the model persists inside the past light cone. 

We thus conclude that the SNIa data does not require a deviation from the Copernican Principle that is non- 
perturbative, at least if we limit ourselves to spherically symmetric perturbations. Indeed, it would be of some 
interest to generalize the present approach to non-spherical perturbations; whether this can be done in practice 
remains to be seen. 



Acknowledgements 

This work was partly supported by the EU 6th Framework Marie Curie Research and Training network "Uni- 
verseNet" (MRTN-CT-2006-035863) and partly by Academy of Finland grant 114419. MM is supported by the 
Graduate School in Particle and Nuclear Physics. MM wishes to thank Teppo Mattsson for useful discussions. 



9 



[1] R. R. Caldwell and A. Stebbins, Phys. Rev. Lett. 100 (2008) 191302 arXiv:0711.3459 [astro-ph]]. 

[2] C. Clarkson, B. Bassett and T. C. Lu, Phys. Rev. Lett. 101 (2008) 011301 arXiv:0712.3457 [astro-ph]]. 

[3] J. P. Uzan, C. Clarkson and G. F. R. Ellis, Phys. Rev. Lett. 100 (2008) 191303 arXiv:0801.0068 [astro -ph]]. 

[4] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517 (1999) 565 arXiv:astro-ph/9812133 

A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116 (1998) 100 9 |arXiv:astro- ph/9805201|. 

P. Astier et al. [The SNLS Collaboration], Astr on. A strophys. 447 (2006) 31 |arXiv:astr~p h/05 10447], 

A. G. Riess et al, Astrophys. J. 659 (2007) 98 [arXiv:astro-ph/0611572] . 
[5] T. Buchert, Gen. Rel. Grav. 32 (2000) 105 l arXiv:gr-qc/9906015| . 
[6] T. Buchert, Gen. Rel. Grav. 40 (200 8) 467 [arXiv :0707.2153l l lgr-qcll. 
[7] S. Rasanen, JCAP 0611 (2006) 003 [arXiv:astro-pri/0607626] . 

[8] S. Rasanen, JCAP 0804 (2008) 026 [arXiv:0801.2692l [astro-ph]]. 

[9] A. Paranjape and T. P. Singh, Phys. Rev. Lett. 101 (2 008) 181101 |arXiv:0806.3497l [astro-ph]]. 

[10] A. Paranjape and T. P. Singh, JCAP 0803 (2008) 023 [arXiv:0801. 15461 [ astro-ph]]. 

[11] E. W. Kolb, V. Marra and S. Matarrese, Phys. Rev. D 78 (2008) 103002 [arXiv:0807.040T1 [astro-ph]]. 

[12] T. Mattsson and M. Ma ttsson, JCAP 0802 (2008) 004 |arXiv:0708.3673l [astro-ph]]. 

[13] T. Mattsson. larXTv:0711.4264l [astro-ph], 

[14] G. Lemaitre, Annales Soc. Sci. Brux. Ser. I Sci. Math. Astron. Phys. A 53 (1933) 51. 

R. C. Tolman, Proc. Nat. Acad. Sci. 20 (1934) 169. 

H. Bondi, Mon. Not. Roy. Astron. Soc. 107 (1947) 410. 
[15] A. Krasinski, "Inhomogeneous Cosmological Models", Cambridge University Press (1997). 

[16] J. Plebanski, A. Krasinski, "An Introduction to General Relativity and Cosmology", Cambridge University Press (2006). 
[17] For a recent review of LTB models, see 

K. Enqvist, Gen. Rel. Grav. 40, 451 (2008) larXiv:0709.2 044l [astro-ph]]. 
[18] M. N. Celerier, Astron. Astrophys. 353 (2000) 63 [arXiv:astro-ph/9907206] . 

H. Iguchi, T. Nakamura and K. i. Nakao, Prog. Th eor. Phys. 108 (2002 ) 809 [arXiv:astro-ph/0112419] . 



D. Garfinkle, Class. Quant. Grav. 23 (2006) 4811 arXiv:gr-qc/0605088 
D. J. H. Chung and A. E. Romano, Phys. Rev. D 74 (200 6) 103507 |arXiv:astro-ph/0608403] 
K. Enqvist and T. Mattsson, JCAP 0702 (200 7) 019 [arXiv:astro-ph/0609120] ~ 
[19] K. Tomita, Astrophys. J. 529 (2000) 26 [arXiv:astro- ph/99052 7'8|7^ 

K. Tomita, Prog. Theor. Phys. 106 (2001) 929 |arXiv:astro-ph/010414r|. 

K. Tomita, Mon. Not. Roy. Astron. Soc. 326 (2001) 287 arXiv:astro-pri/0011484l 



H. Alnes, M. Amarzguioui and O. Gron, Phys. Rev. D 73 (2006) 083519 arXiv:astro-ph/0512006 



H. Alnes and M. Amarzguioui, Phys. Rev. D 74 (2006) 103520 arXiv:astro-ph/0607334 
[20] A. Ishibashi a nd R. M. Wald, C lass. Quant. Grav. 23 (2006) 235 |arXiv:g r-qc/0509108 . 

[21] A. Paranjape, larXiv:0906.3T65l [astro-ph. CO]. 

[22] K. Van Acoleyen, JCAP 0810, 028 (2008) |arXiv:0808.3554l [gr-qc]] . 
[23] M. Sasaki, Mon. Not. Roy. Astron. Soc. 228 (1987) 653. 

[24] E. Barausse, S. Matarrese and A. Riotto, Phys. Rev. D 71 (2005) 063537 [arXiv:astro-ph/0501152]. 
[25] C. Bonvin, R. Durrer and M. A. Gasparini, Phys. Rev. D 73 (2006) 023523 [arXiv:astro-ph/0511183| . 
[26] I. M. H. Etherington, Phil. Mag. ser. 7, 15 (1933) 761. 

[27] G. F. R. Ellis, "Relativistic Cosmology", p. 104 in Proc. School Enrico Fermi, "General Relativity and Cosmology" 

Ed. R. K. Sachs, Academic Press (New York 1971). 
[28] H. Iguchi, T. Nakamura and K. i. Nakao, Prog. Theor. Phys. 108 (2002) 809 i arXiv:astro-ph/0112419] . 
[29] D. J. H. Chung and A. E. Romano, Phys. Rev. D 74 (2006) 103507 |arXiv:astro^p h/0608403|T^ 
[30] C. M. Yoo, T. Kai and K. i. Nakao, Prog. Theor. Phys. 120 (2008) 937 [arXiv:0807.0932l | astro-ph]]. 



[31] N. Mustapha, C. Hellaby and G. F. R. Ellis, Mon. Not. Roy. Astron. Soc. 292 (1997) 817 arXiv:gr-qc/9808079 



